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We analyze the influence of noise for qubits implemented using triple quantum dot spin system. 
We give a detailed description of the physical realization and develop error models for the dominant 
external noise sources. We use a Davies master equation approach to describe their influence on the 
qubit. The triple dot system contains two meaningful realizations of a qubit: we consider a subspace 
and a subsystem of the full Hilbert space to implement the qubit. We test the robustness of these 
two implementations with respect to the qubit stability. When performing the noise analysis, we 
extract the initial time evolution of the qubit using a Nakajima-Zwanzig approach. We find that 
the initial time evolution, which is essential for qubit applications, decouples from the long time 
dynamics of the system. We extract probabilities for the qubit errors of dephasing, relaxation and 
leakage. Using the Davies model to describe the environment simplifies the noise analysis. It allows 
us to construct simple toy models, which closely describe the error probabilities. 



I. INTRODUCTION 

The spin eigenstates of the electron provide one of the 
most natural representations of a qubit - the building 
block of a logical unit in a quantum computer. In recent 
years great progress in the fabrication and the control 
of quantum dots containing only one electron has been 
reported 1 . This progress is essential if a spin-based quan- 
tum computer is to be realized. 

The first proposal of a spin based quantum computer 
used the spin of a single electron in a quantum dot as a 
qubit 2 . In this proposal, single qubit rotations are per- 
formed by pulsed magnetic fields, and a two qubit gate is 
achieved by the Heisenberg coupling of two neighboring 
electrostatically tuned quantum dots. Since electrostatic 
control of a qubit is achievable on much faster time scales 
than control of pulsed external magnetic fields, single 
qubit rotations based on the exchange interaction have 
been proposed 3, . Here an encoded qubit in the Hilbert 
space of two singly occupied spin quantum dots is used. 
The singlet and the spinless triplet level on the two dots 
define the qubit. Manipulation of the singlet-triplet qubit 
has been achieved experimentally 5 ' 6 . 

For universal quantum computation the singlet-triplet 
qubit requires, in addition to intra-dot exchange interac- 
tion, a magnetic field gradient between the two quantum 
dots. It was natural to ask if a different coding of the 
qubit would enable universal computation with the ex- 
change interaction alone; this is realized if the qubit is 
embodied by the states of three singly occupied quan- 
tum dots 7 . The exchange coupling of at least two of the 
three dot pairs should be controllable. Laird et al. 8 and 
Gaudreau et al. 9 have now shown this universal exchange 
control experimentally of the three-electron states in a 
trio of quantum dots. 

The objective of this paper is to explore in detail the 
robustness of this triple quantum dot qubit in contact 
with a realistic set of environments. We have two major 



alternatives to assess, since the spin Hilbert space of the 
triple quantum dot can accommodate a qubit in two fun- 
damentally different ways 10-12 . Recall that three spin- 
1/2 degrees of freedom combine to form four "doublets" 
(total spin 1/2) and four "quadruplets" (total spin 3/2). 
The first approach is to use two out of the four doublet 
energy eigenstates of this system as the qubit levels. To 
manipulate the qubit, we need to control only the sub- 
space spanned by these two states. Consequently this 
qubit is called subspace qubit. But there is a second al- 
ternative: working in the four-dimensional space of states 
with total spin quantum number 1/2, one considers the 
space as a tensor product of two two-dimensional sys- 
tems. One of these two-dimensional systems is taken as 
the coded qubit. This is referred to as subsystem qubit 10 . 
Note that although more abstract, this notion of subsys- 
tem is mathematically identical to that of an ordinary 
subsystem, e.g., the states of one quantum dot in a col- 
lection of many quantum dots. 

A triple quantum dot offers both a subspace and a sub- 
system that are immune against various types of global 
noise 13 ' 14 . Defined in a subsystem a qubit is immune 
against strong collective decoherence, and in a subspace it 
is protected from weak collective decoherence 7 ' 8 . Strong 
collective decoherence is any noise acting globally on the 
triple quantum dot. Weak collective decoherence involves 
just global phase noise. Interaction with the real envi- 
ronment is not simply described by either of these limits 
so our coded qubits will be susceptible to decoherence. 
But the goal is to identify the encoded qubit that is as 
robust as possible against external influences, with the 
longest possible relaxation and dephasing timescales. 

This paper will present calculations of the robustness 
of the subspace and subsystem qubit coupled to realis- 
tic environments for semiconducting spin qubits in triple 
quantum dot systems. We will give a detailed description 
of the qubit implementation and analyze the time evo- 
lution of the noisy qubit. We employ a specific Markov 
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approximation, describing the limit of weak coupling of 
the quantum dot to its surroundings. This model was 
introduced by Davies 15 and will be called Davies Model 
(DM) in the following. We analyze the influence of noise 
and extract error probabilities for relaxation and dephas- 
ing phenomena, as well as for the leakage to other parts 
of the Hilbert space. 

Numerical simulations show that the initial time evo- 
lution can behave quite differently from the long time 
evolution of the qubit. Since we are mainly interested 
in the errors of qubit manipulations that are achieved on 
short time scales, we focus on the description of the initial 
time evolution. We develop an effective master equation 
for the description of the qubit, while removing the in- 
fluence of the environment, using a Nakajima-Zwanzig 
approach 16 ' 17 . The Nakajima-Zwanzig approach espe- 
cially helps to develop a description for the initial time 
evolution. We analyze the initial dynamics in detail and 
describe how error probabilities can be extracted. The 
description in the DM allows us to controllably sort the 
generated dynamics of the quantum dot into groups of 
transition terms. This special structure strongly restricts 
the time evolution of the qubit. Additionally it helps to 
analyze the generated dynamics. We describe the error 
probabilities using a few simple toy models. 

While rather lengthy, we believe that this paper will be 
useful as handbook describing the many possible decoher- 
ence and relaxation scenarios that can arise for the triple 
dot qubit. As more experiments are done to explore the 
various possible encodings of qubits in these systems, the 
results here should serve as a guide to help in arriving at 
the optimal design for making further progress towards 
functioning multi-qubit structures. 

The organization of the paper is as follows: In Sec. II 
we introduce the model analyzed in this paper. We con- 
struct the triple dot Hamiltonian, describe the qubit im- 
plementation and introduce the noise model. In Sec. Ill 
we describe the modeling of real spin qubits. Besides the 
description of the triple quantum dot system, we also 
introduce the noise parameters. In Sec. IV we develop 
the model analysis. We point out how the initial time 
evolution can be analyzed and simplify the description 
method. In Sec. V we analyze the full time evolution of 
the qubit, while in Sec. VI error rates for the initial time 
evolution are extracted. We conclude with a summary 
and an outlook in Sec. VII. 



II. DESCRIPTION OF MODEL 

A. Triple Dot Hamiltonian and Qubit 
Implementation 

The effective Hamiltonian H describing the triple 
quantum dot contains the exchange interaction between 
two neighboring quantum dot pairs. Additionally, an out 
of plane magnetic field is added. Our notations follow 
Laird et al. 8 : 




Figure 1. Layout of the triple quantum dot setup: Each 
quantum dot is occupied with one electron. Neighboring dot 
pairs are tunnel coupled with the coupling strength t. An 
external electric bias is used to occupy either the left (dot 1) 
or the right (dot 3) quantum dot with two electrons. 

U (<ti ■ er 2 - 1) + (er 2 • er 3 - 1) 

»=1,...,3 

J\2 (J23) represents the Heisenberg exchange interac- 
tion between two neighboring dots. It can be changed by 
applying an electric bias on the outer dots. 

The coupling parameters J12 and J23 can be derived 
from a three site Hubbard Hamiltonian describing the 
quantum dot layout in Fig. 1: 

^■Hubbard 

a,s a 

+ t E (°a,«0|8, a + h.C.) (2) 
(a,/3),s 

e a is the single particle energy, U a the coulomb repulsion 
and t the tunnel coupling between neighboring quantum 
dots. For simplicity we take only tunneling into account 
for the left (right) quantum dot with the center dot ((1, 2) 
and (2,3)). Additionally, we assume these tunnel cou- 
plings to be equal. When going from HHubbard to %. 
we take into account single occupation of all three qubits 
((1,1, 1 ^configuration) and use an electric bias to go to 
a double occupied left (right) dot ((2, 0, 1) and ((1, 0, 2)- 
configuration) . For the double occupied states, we con- 
sider only the orbital ground state. 

In analogy to the case of double dots 18 ' 19 , we describe 
all three regimes in an adiabatic basis. The tunnel cou- 
pling causes transitions between the singlet states of all 
charge distributions. In the adiabatic basis we use only 
the ground state of all these charge distributions. We 
eliminate the excited states separately for the charge 
transition to a double occupied left and right dot. One 



3 



arrives at the exchange parameters: 



*- £ V t +l/( £ V t ) > +*" < 4 > 

The bias parameter e lowers the energy of the left quan- 
tum dot for e < 0, while the right quantum dot is fa- 
vored for e > 0. e_ is the bias at which the (1, 1, 1) and 
(2, 0, ^-configurations have the same energy in the ab- 
sence of tunnel coupling (and similarly for e+) 20 . The 
eigenstates of the Hamiltonian % are: 



Qi = 

2 

Qi = 

2 

Ai = 

2 

A'i = 
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where Q = ^/ J 2 2 + J| 3 — J12 J23. For later analysis we 
introduce the notation Wi, indicating the s z = ^ sub- 
space (W € {Q, A, A'}). All remaining eigenstates can 
be obtained by flipping all three spins, to obtain from 
the states Q3 and Wi the corresponding states Q_3 and 
W_i . The eigenenergies of W are: 

E Qk =-k-E z , fce {±5,±Ij (9) 



^A ± , --^(Ji2 + J 2 3-«)T : y (10) 



=-^(Ji2 + J 23 + fi)T^ (11) 

The energy diagram is sketched in Fig. 2. For the 
upcoming analysis we introduce three quantum numbers 
(l,S,s z ), which fully characterize the eigenstates. S de- 
scribes the total spin of the eigenstates. It has the value 
3/2 for all Qfe-states and 1/2 for the remaining ones. 
The s z -quantum number labels the spin projection in z- 
direction. It has values ±3/2 and ±1/2. Furthermore we 
introduce a third formal quantum number I. It distin- 
guishes the Afc (I = 1) and A^ states (1 = 0). 



B. Subspace and Subsystem Qubits 

For our later use we construct a subsystem and a 
subspace qubit inside the 8-dimensional Hilbert space 



spanned by the eigenstates of the triple dot Hamilto- 
nian (1) 11,12 . The basis states of the subspace qubit de- 
fine the computational subspace: span |Ai , A'i |. We 

identify the state Ai as logical "1", A'i as logical "0". 

2 2 
For the subsystem qubit we use the larger subspace 

span \ Ai, A'i, A_i, A'_i > as the computational sub- 

L 2 2 2 2 J 
space. We identify both states Ai and A_i as logical 

"1" (represented by the quantum number 1 = 1). The 
I = states are identified as logical "0". For the sub- 
system qubit the s z -population does not matter. In the 
Nakajima-Zwanzig approach (see Sec. IV A3), we fix this 
population to a constant value. The thermal distribution 
over the Zeeman-split eigenstates will be a reasonable 
choice: 

p s " = e ~^T~ /Tr e (12) 

These two possible ways of defining a qubit are mo- 
tivated by the experimental possibilities for initializing 
the qubit into a controlled state. For the subspace 

qubit, initialization into the A \ -state seems experimen- 

2 

tally achievable. When doing all experiments at high ex- 
ternal magnetic fields (cf. Fig. 2(a)), the s z = — 1 states 
may be avoided; they are far up in energy compared to 
the s z = I states (cf. Eq. (9)-(l 1 ) ) . Initialization into 

A'i can be achieved by coupling two of the three dots 

2 

strongly, creating effectively a strongly coupled double 

quantum dot and an uncoupled single dot. A'i is de- 

2 

scribed by a singlet eigenstate on the double dot, while 
the Ai-state involves triplet eigenstates. Initialization 
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Figure 2. Diagram of the eigenenergies of the exchange 
Hamiltonian of Eq. (1) as a function of the bias parameter 
e with E z > 0. The dashed gray lines are the higher energy 
states that are not treated by the W in Eq. (1). They are 
removed from the 'H.Hubbard by adiabatic elimination, 
(a) shows the case of large external magnetic fields: E z = 
100 /ieV (« 7 T in GaAs). In (b) the external magnetic fields 
are small (E z — 2.5 peV, corresponding to 200 mT in GaAs). 
The dashed orange lines mark regimes analyzed in the analysis 
of Sec. V. 



is now identical to the initialization of the singlet-triplet 
qubit in double dots 19 . The uncoupled single dot needs 
to be in its ground state ( 1 1) ) ■ 

For the subsystem qubit the initialization works out in 
the same way. Here however it does not matter if we ini- 
tialize into the Aj -state or the corresponding A'_ i -state; 

2 2 
both are labeled by I = 0. These states differ only by the 

population of the weakly coupled single dot. For the sub- 
system qubit it is satisfactory to produce a thermal dis- 
tribution between the spin up and spin down states on 
the weakly coupled dot. This property is described by 
the density matrix (12). The strongly coupled quantum 
dot should again be initialized into the singlet state. We 
can accomplish the initialization of the subsystem qubit 



for small and large magnetic fields. 



C. Noise Description 

To study the influence of noise in the triple dot setup 
we first discuss a Lindblad master equation on the 8 di- 
mensional Hilbert space: 

p (t) = (Co + C D ) P (t) = -i [U, P (t)} + C D ( P (t)) (13) 

We will set h — 1 and ks = 1. We add to the coherent 
evolution, given by the Hamiltonian from Eq. (1), a dissi- 
pative Lindblad terms Cd {p (t)) = YIa [A] (p (t)), 
where V [A] (B) = ABA^ - \ (A^AB + BA^A) . An ex- 
ternal bath couples through the operators A to the sys- 
tem. T_4 G R determines the coupling strength. We 
analyze the effects of dephasing and relaxation from ex- 
ternal baths. Dephasing of spin qubits is generated by 
fluctuating magnetic fields parallel to the external mag- 
netic field. Relaxation is generated by fluctuating per- 
pendicular magnetic fields (compare also the description 
in Sec. IIIB). The coupling operators A act either glob- 
ally on the triple dot system or individually on each of 
the dots: 



C g iob (p) = T Z V [Z] (p) + T X V [X] (p) (14) 



£loc (p) 



i=l : 2,3 



(Tf©[4](p) + T^[4](p)) (15) 



With Z = Ei=l,2,3 a l and X = Ei=l,2,3 a x- 

This model represents a specific Markov approximation 
to describe for the time evolution of an open quantum 
system 15 . The procedure for making a Markov approxi- 
mations is mathematically not strict 21 . For our analysis 
we will modify Eq. (14) and (15) by employing a differ- 
ent Markov approximation. Our goal is to make sure that 
the system equilibrates in the long time limit. We will di- 
rectly adopt the description from the paper of Bravyi and 
Haah 22 . They describe a specific Markov approximation, 
originally introduced by Davies for the weak coupling 
limit of a system and a bath 15 . The modified Lindbla- 
dian in the DM is: 



C D (p) = J2h(A,uj)V[A u ] ( P ) 



(16) 



All coupling operators A are decomposed into transition 
terms between equidistant energy eigenstates of the free 
Hamiltonian: 



Here one groups the components of A, so that: 
(A\A U \B) \IE a -E b =uj 



(A\A U \B) - | J [ h gJ wi g e 



(17) 



(18) 
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value 


size 


St/ns 


10 


T K 


10 ^eV (« 125 mK) 


Magnetic field strength 


high: 7 T (E z « 100 ^eV) 
low: 200 mT (E z w 2.5 /xeV) 


exchange interaction parameter 


e+ = |e_| = 500 ^eV 
t = 14 ^eV 


analyzed interval of applied e 


[-500 fj,eV, 500 ^eV] 



Table I. Characteristic values for the analysis of GaAs triple 
quantum dots. The values are used according to the publica- 
tion by Gaudrau et al. 9 . 



B. Transition Rates for Noise Description 

Our model does not contain a microscopic description 
of interactions with the environment. In the DM the 
influence of the surroundings is only modeled by the gen- 
erated transition rates between quantum states (see de- 
scription of Sec. IIC). We will especially focus in our 
analysis on experiments in GaAs. This material has 
been used in all previous experiments on triple quantum 
dots 8 * 9 . 



Hyperfine Interaction 



The rate of quantum jumps h (A, to) G M is set by the 
transition frequencies uj between energy eigenstates of 
the system induced from the bath. It needs to fulfill 
the detailed balance condition connecting positive and 
negative energy differences: 



h (A, -uj) = e T i<h (A, uj) 



(19) 



As shown by Spohn 23 , the Gibbs state is a fixed point of 
the dynamics using the description of the DM: 



Id (e"*) 



(20) 



For simplicity we neglect the tilde on the redefined 
Lindbladian (16) for the remainder of the paper. When 
inspecting the energy diagrams in Fig. 2, one sometimes 
finds that sets of energy levels become equidistant at spe- 
cific exchange interactions, which are not equidistant for 
each e. We won't add these "accidental" degeneracies to 
the construction (16) of the DM. 



III. APPROACH TO MODEL REAL SYSTEMS 

A. System Parameters 

All system parameters to define the Hamiltonian of 
Eq. (1) are matching the triple-dot experiments by 
Gaudreau et al. (compare especially the supplemen- 
tary material 9 ). The parameters are summarized in 
Tab. I. We use a typical experimental time scale St of 
10 nanoseconds 24 . The temperature Tk is set to 125 mK 
(~ 10 fieV). A high magnetic field accounts for the case 
E z 3> Tk, while the low magnetic field case describes the 
opposite limit (E z < T K )- When deriving predictions for 
the initial time evolution, we never go far into the regime 
of double occupied quantum dots. We limit the bias e on 
the interval [e_,e+] of the exchange interaction param- 
eter. The tunnel coupling parameter t is also extracted 
from the paper 25 . 



As in the experiments on single and double quantum 
dots, the nuclear magnetic fields are one major source of 
noise for triple dots. The magnetic moments of the nu- 
clei in GaAs couple through the hyperfine interactions to 
the spin of the electron. Extensive studies of the gener- 
ated dynamics have been carried out in single quantum 
dots 26-29 and double dots 18 . Also very recently a study 
on triple quantum dots appeared 30 . We won't follow the 
arguments of these papers in detail, but will extract tran- 
sition rates for our later analysis. 

In a semiclassical picture the nuclear magnetic mo- 
ments add up to a macroscopic magnetic field 1,31 . Fluc- 
tuations in the nuclear magnetic field are slow compared 
to the precession time of the electron spin in this mag- 
netic field. For the initial time evolution of the electron, 
the nuclear magnetic field can therefore be described as 
static. Due to the large number of spins interacting with 
the electron, one approximates the hyperfine fields by a 
Gaussian distribution with zero mean and a root mean 
square (rms) SE nuc (represented by the distribution func- 



tion / (B) = 



B 

ME*,, 



■-). A typical values in 
0.3 A*eV. 



/2nSE„ nc 

GaAs quantum dots is SE nuc 

We will argue that fluctuating hyperfine fields gener- 
ate locally spin dephasing and relaxation. Its influence 
is reflected by two distinct transition rates, which expo- 
nentially decrease for increasing transition energies uj: 



for to > 
+ 



(21) 



Tk J for oj < 



h {o\ , uj) arises from local magnetic field fluctuations 
at quantum dot i (i = 1,2,3) in the direction of the 
external magnetic field. These operations lead to de- 
phasing of individual spins. h(a l xl uj) describes local 
spin relaxation through in-plane magnetic field fluctu- 
ations (coupling through raising and lowering operators 
<t* = <j\ + (tL 32 )- Quantum jumps between energy levels 
are possible at energy differences smaller than or in the 
range of the hyperfine interaction SE nuc . We will argue 
that both transitions generated from fluctuating out of 
plane magnetic fields and in plane magnetic fields have 
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been observed in previous experiments. We refer espe- 
cially to experiments on single and double quantum dots 
to support the modeling through Eq. (21). 

For single quantum dots a fluctuating mag- 
netic field parallel to the static external field 
leads to dephasing, while a fluctuating perpen- 
dicular magnetic field causes spin flips. In ex- 
periments the time evolution of a single spin 



(Sb„„ c (t) = Tr {<Jp Bnuc (t)) = Tr (<re 



V(0)e 



-iUr 



is measured as the averaged result over many runs of 
the experiment. For each measurement a Hamiltonian 
Ti.B nuc , corresponding to a specific hyperfine field B nuc , 
is determining the time evolution (cf. Fig. 3). The 
measured result reflects the ensemble average over the 
hyperfine field distribution (S (t)) = JS B (t)f (B) dt. 

We use a similar argumentation scheme as in the paper 
by Merkulov et al. 31 . Without external magnetic fields 
the time evolution is completely determined by the fluc- 
tuating magnetic field {Hb — tt &z)- When we calculate 
the time evolution and average it over the magnetic field 
distribution f (B), we see that the component perpen- 
dicular to the magnetic field decreases exponentially: 



(S x (*)) = e 



S x (0) 



(22) 



Since there is no fixed quantization axis of the qubit, 
one expects that all components of the spin decrease. 
The semiclassical analysis predicts therefore a Gaussian 
decay with a time constant (SE nuc ) 1 . The DM describes 
this behavior by a transition rate Y*' x = 8E nuc . This is 
the value of (21) at ui = 0. The energy difference ui 
is determined by the external magnetic field strength. 
In the absence of magnetic fields Tf and Tf must be 
indistinguishable. 

A fixed external magnetic E z and a fluctuating parallel 
magnetic field (% B = fi ', +B cr z ) generates similarly also 
dephasing of the transverse spin components. Now the 
single spin precesses however with the angular frequency 
E z around the z-axis: 



(S x (*)> 



(S x (0) cos (E z t) - S y (0) sin (E z t)) 

(23) 



This is also reflected in the DM. Dephasing transition 
rates for single spins are independent of the energy dif- 
ference. a\ causes only transitions between identical spin 
states, which is reflected by the transition rate h (cr*,0). 
At finite external magnetic fields the relaxations decrease 

with (^jj^ — J ■ We can see this when calculating the 

time evolution for a fluctuating perpendicular magnetic 
field for an initially spin up polarized particle (Zq = 1 



and T-Lb = %cr z 



Ox)- 



Z B (t) = 



E 2 Z + B 2 cos (E z t) 



(24) 



E 2 Z +B 2 

When expanding in (B/E z ) and averaging over the field 



0) 




Figure 3. Semiclassical picture used for the time evolution 
of a single electron spin in a distribution of hyperfine config- 
urations. For every experiment the hyperfine field acts like 
a static magnetic field, which is randomly varying between 
every run of the experiment. 



distribution / (B), we get: 
(Z (t)) « Z 



1 — cos (uit) 



(25) 



Since in the DM relaxation always contain the parame- 
ter h (<7* , uij , we use a Gaussian dependence on the en- 
ergy difference. It describes a quadratic dependence on 



( 8E n uc ) 



for finite uj, when 



(*=)■ 



< 1. 



To extract parameters for h(al,ui\ at finite energy 
differences we can consider double dot experiments. A 
fluctuating local magnetic field parallel to the static ex- 
ternal field causes transitions of the singlet-triplet qubit 
between the singlet state So and the s z = triplet state 
To. On the relevant subspace it acts like a perpendicular 
magnetic field: 



Bal = B(a x ) {SgTo 



(26) 



In the DM the fluctuating parallel magnetic field involves 
quantum jumps between the So and To-states at the en- 
ergy difference of these two levels: h \ a\ , w) . Using the 
same description as for single quantum dots before, the 
transition rates exponentially decrease with the energy 
difference of the levels. 

Transitions between the s z = 0-states with the s z ^ 
states are possible through local raising and lowering op- 
erators. They are generated through perpendicular mag- 
netic field fluctuations h (a l x , uA at the energy difference 
uj of the s z = and s z ^ states. The transition rates 
are highly sensitive to uj. This result was also observed in 
experiments on double dot qubits. The transition rates 
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Figure 4. Hyperfine interaction causes dephasing and relax- 
ation through local magnetic field fluctuations on each quan- 
tum dot. Dephasing and relaxation is generated by the tran- 
sition rates h \a\j '<j % x ,iJ) according to Eq. (21). The marked 
points refer to results from two qubit experiments, described 
in the main text. 



of (20 ns) at zero magnetic fields and (f 50 ns) at 
100 mT in our model (cf. Fig. 4) match approximately 
these results 6 . For weakly coupled singlet-triplet qubits 
the energy difference w is determined directly by the ex- 
ternal magnetic field. 

Interestingly we automatically describe in the DM the 
large transition rates at the crossings of levels with dif- 
ferent spin quantum numbers. For double quantum dots 
the Sq-T + crossing is extensively used in experiments 5 . 
In triple-dots we will observe doublet-quadruplet cross- 
ings, which were analyzed in the experiment by Gaudreau 
et al. 9 . In experiments a large enhancement of transition 
rates near these level crossings were observed, with tran- 
sition rates in the nanosecond regime. These transition 
rates decrease quickly when going away from the level 
crossing. 



2. Interactions with Phonons 

Additionally, relaxations are generated by interactions 
of the quantum dot with electric fields, e.g. from polar 
phonons. However, direct transitions are forbidden be- 
tween different spin states of the same orbital level by 
the dipole selection rule. They must be mediated by an- 
other process like spin-orbit interaction. In single quan- 
tum dot experiments, relaxation times of about 1 s at 
magnetic fields of 1 T and 0.5 ms for a magnetic field of 
5 T have been identified 33 . The scaling law of uj 3 E% gov- 
erns transitions by piezoelectric phonons between Zee- 
man split quantum dot eigenstates for single qubit ex- 
periments. The phonon energy must match the energy 
difference of the states, resulting in the w 3 -scaling law. 
The Zccman split eigenstates are mixed by spin orbit in- 
teraction, which causes the E% dependence. Experiments 




Figure 5. Relaxation rates h (er^.,u;) generated by piezoelec- 
tric phonons according to Eq. (27). Experimentally observed 
time scales from qubit experiments 33 are included in this pic- 
ture. 



mechanism constant 


value 


local dephasing and relaxation y Z ,x 
through hyperfine fields 


i 

20 ns 


local spin relaxation „ x 
through phonons 


2 • 10" 6 * 

s pev J 



Table II. Transition mechanisms describing noise on the triple 
quantum dot and estimates of parameters. The different 
mechanisms are described in the text. 



and theory suggest therefore relaxation rates modeled as: 



h(a l x ,uj) = 



1 



(27) 



This should also apply for triple dot setups. The coupling 
operators are again local spin flip operators a x = a'T+crT' 
at quantum dot i. A picture of the generated transition 
rates, including the results from single qubit experiments, 
is shown in Fig. 5. 

For double dot experiments in the weak coupling 
regime, transition rates have been identified that are con- 
sistent with this picture 34 . At high bias (readout regime) 
other effects cause transitions between different charge 
states 19 . We do not include these effects in our model 
since we do not analyze the readout regime (e > e + or 



3. Summary and Description of Noise Regimes 

In the previous analysis we identified two influences 
of the qubit environment, which should be most impor- 
tant for triple dot experiments. Hyperfine interactions 
cause large error rates through the random distribution 
of their magnetic moments. Direct spin flips are gener- 
ated by phonons, which couple to the spin indirectly via 
spin orbit interaction. We want to point out that charge 
fluctuations might be an additional noise term in triple 
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Local phase noise h {a\,uj) generated from 
Regime 1 fluctuating hyperfine fields. Strong influence 
at small energy differences. 

Local spin relaxation h (0^,0;) from the interaction 
Regime 2 w jth phonons. Strong influence at large energy 

differences and large external magnetic fields. 

Local spin relaxation h (cr*,u;) from the interaction 
Regime 3 with hyperfine fields. Dominant relaxation 

mechanism at small energy differences, independent 

of the external magnetic field. 

Table III. Dominant noise regimes for qubit experiments. 

quantum dot experiments. Through the exchange inter- 
action electrostatic interactions are used to manipulate 
the qubit. Therefore charge fluctuations will be impor- 
tant. This problem was already pointed out for double 
quantum dots 35 . Especially for strong electrostatic bias, 
charge dipoles are created and charge fluctuations will 
gain influence. Here we study only the weak bias regime 
and do not take into account charge fluctuations. 

When we analyze the qubit, we will find that the de- 
scribed noise mechanism will be relevant for distinguish- 
able parameter sets. A summary of the three most impor- 
tant regimes is given in Tab. III. For phase noise we de- 
scribe only the influence of hyperfine interaction through 
the transition rates h {a\, w) from Eq. (21). This interac- 
tion will cause major errors for energy differences smaller 
and in the range of the rms of the hyperfine energy (cf. 
Fig. 4). This regime will be the most critical one for ex- 
periments in semiconducting spin qubits. We will call it 
"Regime 1" in the following. 

Local spin relaxation can cause large transition rates in 
two completely different regimes (see transition rates in 
Fig. 6). For large energy differences and external mag- 
netic fields the interaction with phonons will be domi- 
nant, while at small energy differences the hyperfine in- 
teraction will determine the relaxation process. Since 
these two effects are important in different parameter 
ranges, we can easily separate their influence. We will 
call the two parameter ranges "Regime 2" and "Regime 
3" respectively. 

IV. PREPARATION FOR MODEL ANALYSIS 
A. Description of Initial Time Evolution 

1. General 

In qubit experiments one is usually interested in the 
time evolution of the qubit on short time scales. With 
the Nakajima-Zwanzig approach one can construct an ef- 
fective master equation for the initial time evolution of 
the qubit 16,17 . The "common" master equation describes 
the time evolution of the full system (with its multiqubit 




0) [pteV] 



Figure 6. Density plot of the transition rates for local spin 
relaxation h (o 1 x ,lu) as a function of the energy difference to 
and the external magnetic field E z . The influence of hyper- 
fine interaction (Regime 2) and phonon interaction (Regime 
3) can be separated due to the different parameter ranges 
when the effects are dominant. The hyperfine mechanism is 
independent of E z . 

Hilbert space) in a first order differential equation. Us- 
ing the Nakajima-Zwanzig approach one can reduce this 
equation to the relevant part of the Hilbert space describ- 
ing just the qubit. In general the problem of solving this 
lower dimensional equation is not simpler than solving 
for the dynamics of the full system. But with a few addi- 
tional assumptions we can simplify this lower dimensional 
equation. 

As a first step one identifies a relevant part of the 
Hilbert space H rel C H, which is used to define the qubit. 
One defines a linear map V, which constructs from the 
full density matrix only the relevant part: 

p rel (t)=Vp(t) (28) 

We need only two properties for the map V to be phys- 
ically meaningful. First of all the map should act on 
the relevant part of the density matrix like the identity 
operation. One disposes the condition: 

V 2 =V (29) 

Secondly an observable T on the relevant part of the 
Hilbert space should be described in the the same way 
by the reduced density matrix V p (t) and the full den- 
sity matrix p(t). We obtain this physical property by 
requiring: 

Tr ( TV*) = Tr (30) 

• represents an arbitrary element of Liouville space. Fi- 
nally for our later purpose we also add a third charac- 
teristic. Initially the qubit is decoupled from the sur- 
roundings, which gives Vp(0) = p(0). This requirement 
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is equivalent to the criterion to initialize the qubit into a 
controlled state. 

With these three assumptions we will rewrite our Lind- 
blad master Eq. (13): 

p(t)=Cp(t) (31) 

For the upcoming analysis C can consist of a coherent 
time evolution £ (•) = — i \Hi •] an d it may also include 
a dissipative Lindblad term Cd (•) = ^2 a [A] (•). 
One can exactly rewrite the master Eq. (31) for the rel- 
evant part p rel (t) with a time retarded equation 16,17 : 

p rel (t) =T£Tp rel (t) (32) 

+ / dt'VCQe^^-^QCVp™ 1 (t') 
Jo 

Eq. (32) is called the Nakajima-Zwanzig equation. We 
have introduced the projector Q = 1 — V. 

To describe the initial time evolution, we divide the 
full Hilbert space into a relevant part A and an irrele- 
vant part B. Cq = C A + C B describes the time evolution 
of A and B individually. C\ connects A and B. The time 
evolution should be dominated by Co, while C\ describes 
only a "small" term. In second order Born approxima- 
tion, we keep terms containing C\ up to second order. 
The Nakajima-Zwanzig equation in second order Born 
approximation reads: 



superoperators. They describe transitions between two 
of these Liouville subspaces. The superoperator C p has 
only components connecting superstates from P and P: 

C P = \P)L PP (P\ (35) 

Cq never acts on the relevant subspace: 

C Q = \Q)L QQ (Q\ (36) 

The remaining superoperator Cpq does not just couple 
the subspaces P and Q. It also has contributions to the 
off-diagonal terms of the density matrix: 

C PQ = \A)L AB (B\ (37) 

A,B€{P,Q, + ,-} 
AB${PP,QQ} 

Using this notation one can rewrite Eq. (33) for the 
linear map Vp from Eq. (34). We arrive at a master 
equation on the relevant subspace P: 

p p (t) =L P pp p (t) + f dt?Lp Q e Le >°( t - t ')L Q pp p (t) 
Jo 

+ (L P+ L +P + Lp_L_p) [ dt'p p {t') (38) 

Jo 



p p {t)={VC A V + VC AB V)p p {t) (33) 

+ f dt'V£ AB Qe QiCA+CB)Q ( t - t ')Q£ABVp p (0 
Jo 

2. Subspace Qubit 

To define the subspace qubit one needs to project out 
parts of the Hilbert space, i.e. one uses a map Vp made 
up of projectors. The parts A and B from Sec. IV A 1 are 
subspaces of the full Hilbert space. We we will call them 
P and Q respectively. Vp constructs from p (t) € H ~ C d 
the relevant density matrix on subspace P. p p (i) is only 
nonzero in H p ~ C 2 (2 < d). The linear map Vp can 
be constructed to keep from full density matrix only the 
relevant components: 



V P : p = 




Qp is implicitly defined as 1 — Vp. To rewrite Eq. (33) 
for V = Vp and Q = Qp, a bra-ket notation of superop- 
erators turns out to be very useful (for an introduction 
see the book of Blum 36 ). We use round brackets for su- 
perstates in Liouville space. The superprojectors \i) (i\ 
project onto the corresponding part of the density ma- 
trix. They also divide Liouville space into four subspaces, 
which we label by i, i <E {P, Q,+,— }. We can rewrite 
all superoperators in this notation and identify projected 



3. Subsystem Qubit 

In general one is interested not only in dividing the 
Hilbert space into two subsystems, but also in defining a 
subsystem inside a subspace of the full Hilbert space: 



H= \U S ®H B ©H Q (39) 
V Hp / 

We need this approach for the definition of the subsys- 
tem qubit (cf. Sec. IV A 2). Here we first project on a four 

dimensional subspace P = span \ A i , A'i , A_i , A'_i \. 

L 2 2 2 2 J 

Inside the subspace P, we identify a two dimensional sub- 
system S to define the qubit. For the subsystem qubit, 
the subsystem S is specified by the formal quantum num- 
ber I. The irrelevant subsystem B is characterized by the 
s z -quantum number (cf. Sec. IV A 2). 

The projection of the master equation on the P- 
subspace works in the same way as described in 
Sec. IV A 2. We only need to use a projector V = Vp 
on a four dimensional subspace. We now study the mod- 
ification of the effective master equation due to the intro- 
duction of the subsystem S in the F-subspace. We start 
with the master equation p(t) = C (t) p (t) with a time 
dependent superoperator defined in Eq. (38): 

£(t) = L P p + T(t) (40) 

with: 



10 



T(t) P p (t) = 



= / dt' 
Jo 



LpQe 



>(*-*')l 



Lp+L+p + Lp-L-p I p p (i') 



(41) 



S(t-t') 



T (t) integrates the density matrix over all past times. 
To describe the evolution on the subsystem, one uses a 
linear map V = Vs consisting of a partial trace: 

Vs:p P {t)^p*Tr B (p p {t)) (42) 

p s (t) 

The linear map Vs fulfills especially the properties (29) 
and (30). It extracts from the density matrix of the sub- 
space p p (t), the density matrix of the subsystem p s (t). 
It should be emphasized that we exclude entanglement 
between the systems S and B through the choice of the 
map in Eq. (42). We fix the subsystem B to a static value 
Pq. The effective master equation for the subsystem S 
can be rewritten for the time dependent superoperator 
(40), as shown by Fick and Sauermann 16 ' 17 : 

±pS(t)=r s £(t)V s p s (t) (43) 

+ /* dt'VsC (t) Q S V (t, t') Q S C (t) V s p s (f) 
Jo 

with: 

±V(t,t') = Q s C(t)V(t,t') (44) 

In the analysis of triple quantum dots we try to extract 
errors for the initial time evolution. For this purpose, we 
can rewrite the effective master equation (43) and (44) 
for the description of short times. We divide the Lind- 
blad operator Lpp from Eq. (40) into a part which acts 
just on the qubit subsystem S (Cs) or the irrelevant sub- 
system B (Cp) individually. The remaining dissipative 
term is identified by the operator Csb- Csb should be 
small compared to the Cs and Cb- In second order Born 
approximation we get the effective master equation: 

j t P S (t) =V S (Cs + C SB +T (*)) Tsp s (t) (45) 

+ f dt'VsCsBe^+^^QsCsBVsP 3 (f) 
Jo 

B. Simplification of the Analysis 

1. Rotating Frame 

When analyzing the master equation, we are interested 
in the deviation of the qubit evolution from the free evo- 



lution Co due to the dissipative Lindblad term Cp (cf. 
noise description in Sec. II C). It is therefore meaningful 
to go for the further analysis to a rotating frame with 
respect to the free Hamiltonian: 

p (t) -+ p rot (t) = U rot (t) p {t) U\ at (t) (46) 

with U rot (t) — e mt . This leads automatically to a redef- 
inition of the Lindbladian: 

Co + Cp> > C rot (47) 

Due to the general conditions of the DM (16), the Lind- 
bladian in the rotating frame equals the original dissipa- 
tive Lindbladian: C rot = Cp,. Cp, consists of a sum of 
terms, in each of which the coupling operators A u appear 
twice - once as hermitian conjugate. When the coupling 
operators are written in the eigenenergy represen- 
tation of the free Hamiltonian, to each entry a complex 
argument e luJt is added. The phase ui represents the en- 
ergy difference of the states that the coupling operators 
A u connect. Since all Lindblad operators are grouped 
to couple only equidistant energy levels, these complex 
factors cancel out. 

2. How to Extract Error Rates'? 

To describe the initial time evolution we will see that 
it is sufficient to use a set of 7 just parameters (see also 
the description in Appendix B). This initial evolution is 
a trajectory on the Bloch sphere with full rotation sym- 
metry around the z-axis and reflection symmetry to any 
plane containing the z-axis. Consequently we can restrict 
all our analysis to one plane (e.g. x-z-plane). Addition- 
ally, the trajectory starting on one the north or the south 
pole of the Bloch sphere will be strictly restricted to the 
z-axis. 

These symmetries are very specific to the analysis of 
the problem in the DM (16) and the specific form of the 
quantum jump terms (see Eq. (14) and (15)). The dis- 
sipative terms of the DM are collected to generate tran- 
sitions between equidistant energy levels through the su- 
peroperators V [AJ[. This picture will prove to be very 
helpful to explain the symmetry of the trajectory. The 
Lindblad operator is the generator of time evolution for 
the density matrix. It maps the initial density matrix to 
the density matrix at some later time: 

p(St) = e CSt p(0) (48) 
= (l + CSt + C 25 ^- + ...^j p(0) (49) 
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Eq. (49) makes it clear that all possible combinations 
of superoperators V [A^] V [B^] V [C„] . . . will act on the 
initial density matrix to generate the density matrix at 
some later time. For the subspace qubit we start with a 
density matrix: 

( Og + Zg Xg-iYg \ 
Xg+iYg Oq-Zq 1 
2 2 / S=\,s z = 

e(o 2 ) s4A= _.e(o 4 ) s=i 



Initially there is no population in the subspace spanned 

by ! | and in the quadruplet subspace. It can 

be proven easily (see below) that the action of any com- 
bination of quantum jumps on the density matrix will 
lead to a density matrix of this form: 



(50) 



C n (p(0)) 



axOo + a 2 Z Q a 3 (X - iY ) \ i 
a 3 (X + iY ) a 4 O Q + a 5 Z I \ 

( I1O0+I2Z0 

73O0 + 74^0 

75O0 + 76^o 

V 



fcOo + foZo (3 3 (X -iY ) 
/3 3 (X +iY ) /3 4 O a + /3 5 Z 



\ 



(51) 



77O0 + 78^0 / 



The coefficients 7^) are real numbers represent- 

ing the action of quantum jumps between energy levels. 
Inspecting the density matrix of Eq. (51), it is clear that 
the projected part on the qubit subspace will have the 
same ratio between the x- and y-polarization as the initial 
density matrix. The trace evolution and the z-evolution 
is however dependent on the initial z-polarization of the 
qubit. Since this finding is true for all summands of 
Eq. (49), it is also true for p(St). Given these restric- 
tions on the generated density matrix, the trajectory on 
the Bloch sphere will have the specific form described 
earlier. 

We point out how to prove Eq. (51) with some easy cal- 
culations. All quantum jump transitions can be grouped 
into two sets. First there are transitions involving only 
the computational subspace. They can represent pure 
relaxation (model 1 in Appendix C 1) or pure dephasing 
(model 2 in Appendix C 2) . These models generate tran- 
sitions in the computational subspace via the coupling- 
operators (cf. Eq. (17)): 




(52) 



with real coefficients /3j. An initial density matrix po = 
( ai O + a 2 Z a a 3 (X zF ) \ ^ w gtructuraUy 
y a 3 (X + iY ) a 4 O +a 5 Z J 

the same form after the action of one dissipative term 



(T> [A] (po))- Only the constants on will be modified. 

Transitions involving the remaining Hilbert space will 
have again two distinct features. First there are quan- 
tum jump terms involving just transitions between two 
energy levels. One can calculate the action in a three 
dimensional Hilbert space spanned by the two qubit lev- 
els and the coupled external level. An easy calculation 
shows that the structure of an initial density matrix: 

^ a 1 OQ + a 2 Z a 3 (X — iY Q ) ^ 

pa = a 3 (X + iY ) a 4 Q + a 5 Z 

V a 6 O + a 7 Z J 

(53) 

will remain unchanged. Again only the constants on will 
be modified to different real numbers. 
Secondly there are quantum jumps involving more than 
3 energy levels. They can be made up of all the tran- 
sitions introduced so far, but acting on separate sub- 
spaces. Consequently they also preserve the structure 
above. Otherwise they couple the computational sub- 



space to the (S 1 



2' ' 



2) -subspace involving cor- 



related quantum jumps between the same Z-eigenstates. 
These transitions preserve also the off-diagonal elements 
of the density matrix on the computational subspace. 
They never mix diagonal and off-diagonal elements. 

The same result is obtained for the subsystem qubit. 
In fact all arguments will be identical, since the initial 
density matrix of the subsystem qubit is already in the 
form of Eq. (51): 
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P(0) 




e (o 4 )c 



3. Symmetry of Phase Noise 

We want to point out a key symmetry for phase noise, 
which will simplify the upcoming considerations. The 
action of phase noise through the coupling operators a\ 
(i = 1,2,3) will have an equal effect on the s z = \ 
and the s z = — \ subspace (involving also the quadru- 
plet levels). It will just mix within these subspaces but 
never couples subspaces of different s z -quantum number. 
Furthermore the corresponding matrix elements in the 
s z = \ and the s z = — | subspace are, up to a sign, iden- 
tical. This can be understood by the symmetry opera- 
tion which flips the spins on all dots Ufu p . It transforms 
a state from the s z = +| subspace to the corresponding 
s z = — | subspace and vice versa. It also adds a sign to 
a\. This proves that: 

(W 1/2 \ai\V 1/2 ) = - (W_ 1/2 \ai\V_ 1/2 ) (54) 

for W, V G {A,A',Q}. Since in every dissipative term 
these matrix elements appear twice, the factor " — 1" 
drops out. This symmetry was also identified in the pa- 
per by Ladd, which however did not connect it to the 
underlying symmetry operator 30 . 



4- High Symmetry Regimes 

We point out high symmetry regimes of the qubit 
Hamiltonian, which help us to understand the results 
later on. First of all without bias we will have effec- 
tively a spin zero or spin one particle from the electrons 
of the outer two dots coupled to a spin- 1/2 particle on the 
middle dot. This can be seen easily when noticing that 
without bias J\ 2 = J23 = J- The exchange interaction 
part simplifies to: 

3\ 2 <J\ ■ cr 2 + J 2 3(T 2 ■ er 3 = J cr 2 ■ (ai + er 3 ) (55) 

We can construct eigenstates of the triple-dot Hamilto- 
nian (1) from spin-1/2 eigenstates on the middle dot and 
singlet-triplet levels on the outer two dots. 

Secondly in the case of large positive (negative) de- 
tuning the exchange interaction parameter J 2 3 (^12) is 
dominant. We can ignore the coupling of one dot. Hence 
the model describes a strongly coupled double dot and 



an uncoupled spin-1/2 level. The eigenstates for the ex- 
change interaction Hamiltonian (1) can again be con- 
structed from the singlet-triplet eigenstates of the double 
dot and the single electron eigenstates of the uncoupled 
single dot. We summarize all eigenstates in Fig. 7. The 
corresponding s z = —^-states can be obtained by flip- 
ping all spins. These eigenstates are agreeing with the 
limits of Eq. (6)- (8). 

Additionally, we can understand the action of lo- 
cal noise more easily, when additional symmetries are 
present. First of all noise on dot 1 is equivalent to noise 
on dot 3 when changing the sign of e. This property is 
only true, because in our analysis the tunnel coupling be- 
tween the dot pairs 1 and 2 is identical to the coupling 
between dots 2 and 3. Additionally, we use |e + | = |e_|. 
We therefore never analyze noise on dot 3 individually. 
For the same reason local noise on dot 2 is equivalent for 
positive and negative bias. Secondly, for large negative 
detuning it does not matter if the noise is acting on the 
first or the second dot. This result is just a consequence 
of the situation described earlier. For large negative de- 
tuning we couple dot 1 and 2 strongly, while the third 
quantum dot is effectively decoupled. 

V. FULL TIME EVOLUTION 

Before describing the initial time evolution in detail we 
will go one step back and analyze the full time evolution 
of the qubit under the influence of noise. We will then 
be able to show that the Nakajima-Zwanzig equation de- 
rived in Sec. IV A describes the qubit evolution well on 
the time interval [0, <5£]. 

We numerically calculate the time evolution for a sub- 
space and subsystem qubit on the full 8-dimensional 
Hilbert space. We solve the full master equation in the 
rotating frame: 

p(t)=£ rot p(t) (56) 

The Lindbladian C rot is given by £ D of the DM from 
Eq. (16), as described in Sec. IV B 1. From the full 
density matrix we calculate the time evolution of the 
qubit's population O (t) = Tr (Vp (t)). Additionally, 
we describe the trajectory on the Bloch sphere P (t) = 
(X (t) ,Y(t),Z (t)), with P t (t) = Tr (a.Vp (t)). For the 
subspace qubit we use the map Vp from Eq. (34). It 
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Figure 7. High symmetry regimes of the exchange interaction 
Hamiltonian (1) in the limits of high bias (e = ±00) and 
no bias. The s z = | eigenstates describe always composite 
system of two spin-1/2 levels coupled to one spin-1/2 level. 



projects from the full 8-dimensional Hilbert space on 
the 2-dimensional Hilbert space describing the qubit. 
To define the subsystem qubit, we use the combina- 
tion V — VsVp. Vp is a projective map, similar to 
the construction of the subspace qubit. For the sub- 
system qubit it projects however on the 4-dimensional 

subspace span Ai, A'i , A_ 1 , A'_ k > (see description in 
L 2 2 2 2 J 

Sec. IV A 3). Vs is defined in Eq. (42). We compare 
the results with the solution of the Nakajima-Zwanzig 
equation in Born approximation from Sec. IV A. For the 
subspace qubit we use Eq. (38), for the subsystem qubit 
Eq. (45). We solve both Nakajima-Zwanzig equations 
numerically. We keep noise terms in first (second) order 
for the first (second) order Born approximation. 

We discuss in the following the subspace and sub- 
system qubit separately. Even though we discuss the 
time evolution quite generally, it will turn out that 
the environment influences the qubit evolution domi- 
nantly through transition rates from three different noise 
regimes. We have specified them already in Sec. IIIB3 
and refer to them in the following (cf. Fig. 6). 



A. Subspace Qubit 



The subspace qubit is defined on the subspace 

span |Ai , Aj_ j (see description in Sec. II A). We discuss 

it using all parameters from Sec. Ill A. It constitutes a 
decoherence free subspace with respect to global phase 
noise T> [Z] (weak collective decoherence, cf. Sec. I). 
Other global noise (through the action of the Lindblad 
operators V [X] and V[Y], with X 7 Y = £i=i,2,3 <4, y ) 
will however lead to leakage from the computational sub- 



space. This leakage will separately bring the I = 1 and 
I = states to thermal equilibrium. But since the sub- 
space qubit can only be operated in the limit of large 
external magnetic fields, the resulting leakage will be neg- 
ligible. For E z = 100 /^eV and Tk = 10 fieV only a frac- 
tion e~"fe / (l + e~^^j w 4.5 • 10~ 5 of the probability 

will leak out of the computational subspace. In any case 
we argued in Sec. IIIB that the external environment 
couples dominantly through local interactions. 



1. Regime 1 

Local phase noise turns out to be critical for the sub- 
space qubit. We discuss in the following an unbiased 
triple dot (e = 0). Additionally, we will include phase 
noise generated from fluctuating hyperfine fields (cf. 
h (<7*,w) from Eq. 21). We will include phase noise on 
all three quantum dots. The parameters Tf = (20 ns _1 ), 



30 ns- 



and T3 



(40 ns I describe the com- 



mon noise strength for fluctuating hyperfine fields (see 
Tab. II). Local phase noise generates large transition 
rates only at small energy differences (Regime 1, cf. 
Tab. III). 

A qubit initialized in the s z = ^-subspace will only 
undergo transitions within this subspace. Since the en- 
ergy differences of all three energy levels are comparable 
to the thermal energy (see Eq. (9)-(ll)), a considerable 
part of the population leaks into the Qi -subspace. The 
finite energy difference of the two qubit levels leads in 
the long time limit to a finite polarization on the z-axis: 

•Zoo- 

We show in Fig. 8 the time evolution for the unbiased 
subspace qubit (e = 0). It is initially in a pure excited 
state (P (0) = (0,0,1), see Fig. 8(a)) or in a superposi- 
tion of Ai and A'i (see Fig. 8(b)). For the initial time 

2 2 

evolution (orange interval) the Born approximation in 
Eq. (38) gives a highly accurate description. One can 
see that the thermalization happens on a microsecond 
timescale. Simple arguments can be used to predict the 
long time limit of the time evolution. The arguments in 
Appendix A indicate Ooo rs 0.66 and P m « (0, 0, —0.03). 
The small value of the polarization in the long time limit 
(Zoo —0.03) reflects that the thermal energy Tk is 
large compared to the energy splitting of the qubit levels 
(Tk = 10 ^eV; cf. energy diagram in Fig. 2(a)). 



2. 



nme 2 



Local spin relaxation of the subspace qubit at high 
magnetic fields will predominantly depopulate the qubit 

by the transition from the qubit eigenstates Ai and A'i 

2 2 

to states with different s^-quantum numbers. We simu- 
late the time evolution for the unbiased subspace qubit 
(e = 0) at E z — 100 ^eV (see Fig. 9). Large transi- 
tion rates are generated from the interaction with piezo- 



14 



£0.8 
O 



0.6 
1 



S 0.8 

N 

0.6 



10 



20 



-Lindblad equation 
■ first order Born 
■second order Born 
• t — > oo 



30 



40 




50 



1 

0.5 

^X^V. 


t[/j,s] 


1 2 





10 20 30 40 
t[ns] 

(a) P(0) = (0,0,1) 



50 



0.5 



50 



100 



150 



200 



250 





— Tjinnhlfiri pmiflti nn 

1 ' 1 LL L.1 \-) LCAiKA. v3 V_J LI CL Li 1 \J ±1 

■ ■ first order Born 
■■■second order Born 

— t — > oo 


.1 . 1 






50 



100 150 200 250 

t[fj,s] 

(a) P(0) = (0,0,1) 




20 30 
t[ns] 



0.5 



50 



100 



150 



200 



250 



N -4 





— Lindblad equation 

- - first order Born 
■■■second order Born 

- - 1 — > oo 






250 



(b) P(0) = (0,1,0) 



(b) P(0) = (0,1,0) 



Figure 8. Time evolution of the subspace qubit without elec- 
tric bias (e = ^ieV) at E z — 100 £teV. The qubit is subjected 
to local phase noise generated by fluctuating hyperfine fiefds 
(Tf = (20 ns- 1 ), XI = (30 ns- 1 ), T z s = (40 ns' 1 )). The 
orange region marks a typical time interval for qubit exper- 
iments: [0, 10] ns. Blue lines represent the results from the 
first order Born approximation, while red lines are calculated 
in second order Born approximation. The insets show the 
evolution on a longer time scale. Green lines represent the 
long time limit of the evolution as discussed in the text. 



Figure 9. Time evolution from local spin relaxation for the 
unbiased subspace qubit generated by the interaction with 
piezoelectric phonons for e = and E z = 100 fieV. The 
transition parameters are chosen to ~E X =2-10 

1.5 • lO" 6 ^^ and Hf = 1 ■ lO" 6 ^^. 
solution of the full master equation (black lines), the first 
(blue lines) and second (red lines) order Born approximation 
are identical. The qubit depopulates completely in the long 
time limit. 



6 1 
s ^eV 5 ' " 

The numerical 



electric phonons. These rates are highly enhanced for 
high energy differences and large external magnetic fields 
(Regime 2. cf. Tab. III). We use the following transition 
parameters: 3« = 2 • lO" 6 ^, S* = 1.5 • 10~ 6 ^ 

and 3| = 1 • lO- 6 ^^ (cf. Tab. II). At low electric 
bias we will only see phonon generated transitions to the 
Qa-state. 

2 

Local spin relaxation from fluctuating hyperfine fields 
can be neglected. This mechanism is only dominant 



for small energy differences (Regime 3, cf. Tab. III). 
For large external magnetic fields, the unbiased subspace 
qubit does not have a state of different S2-quantum num- 
ber, which is close to the qubit states (cf. Fig. 2). 

It can be seen that the transitions empty the qubit 's 
population completely, but the time evolution is very 
slow. Overall, phonons generate considerable changes 
only on microsecond timescales. The first order Born 
approximation is already sufficient for the description of 
the qubit dynamics. 
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B. Subsystem Qubit 

A subsystem qubit is implemented in the formally 
introduced Z-quantum number. We define the A±i/2 
and A^_jy 2 -states as the qubit levels (see description in 
Sec. 11 A). It operates in a decoherence free subsystem 
to all interactions acting globally on its defining subsys- 
tem (strong collective decoherence, cf. Sec. I). All global 
noise mechanisms will hence be irrelevant. In any case, 
we argued in Sec. II C that local interactions dominate 
the noise properties in triple quantum dot experiments. 

Local phase noise will have the same effect on the sub- 
system qubit and the subspace qubit. As described in the 
symmetry discussion in Sec. IV B 3, phase noise will act 
separately on the s z = ±| -subspaces. It will however 
never mix them. Additionally, the action is the same 
on both subspaces. Since the transition rates from the 
interactions with nuclear spins are assumed to be inde- 
pendent of the magnetic field strength (cf. Eq. (21)), the 
time evolution is also identical. Large transition rates 
are generated for energy levels that are close in energy 
(Regime 1, cf. Tab. III). 

For local spin relaxations the description of the subsys- 
tem qubit is comparable to the subspace qubit for large 
external magnetic fields. We can especially see, that the 
s z -distribution (12) is close to a pure s z = +|-state. 



1. Regime 3 

Local spin relaxation will gain importance for small 
magnetic fields. In the low bias regime (i.e. small 
|e|) relaxation rates generated from the interaction with 
phonons are small. We can see this by inspecting the 
transition rates h (<rj. , w) for small magnetic fields E z 
(cf. Eq. (27)). On the other hand, fluctuating hyper- 
fine fields can strongly mix different s z -states. Especially 
at the points of level crossings this effect will be critical. 
The relaxation effects through nuclear magnetic fields are 
highly enhanced at small energy differences (see transi- 
tion rates h(a x ,u>) in Eq. (21)). This dominant noise 
mechanism is summarized in Regime 3 of Tab. III. 

We simulate the qubit evolution at e = 354.6 neV and 
E z = 2.5 ueV. Here two different doublet levels and also 
doublet and quadruplet levels are close in energy (see 
orange line in Fig. 2(b)). We simulate the time evolu- 
tion in the rotating frame (cf. Eq. (56) with C rot = Co 
from Eq. (16)) and compare it with the results from 
the Born approximation in Eq. (45). The results are 
shown in Fig. 10. We use following transition parame- 
ters to model the hyperfine interaction: Tf = (20 ns _1 ), 
Tf = (30 ns- 1 ), Tf = (40 ns" 1 ) (cf. Tab. II). 

Local spin relaxation, generated from fluctuating hy- 
perfine fields, mixes especially two subspaces. First 

of all, the subspace span |<3i/2, A1/2, A'_jy 2 | is mixed 
strongly. Secondly, also transitions between the levels 
A'jy 2 and Q3/2 are strong. In the long time dynamics, 
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Figure 10. Time evolution of the subsystem qubit at 
E z — 2.5 fieV through local spin relaxations at e = 354.6 fieV 
(compare energy diagram in Fig. 2(b)). The major influence 
is from local spin relaxation, generated from hyperfine in- 
teractions. We use Tf = (20 ns -1 ), Tf = (30 ns -1 ) and 
Tf = (40 ns" 1 ). The orange region marks a typical time scale 
of qubit experiments. Here the second order Born approxi- 
mation and the results from the full master equation match 
very closely. The insets show the long time evolution of the 
qubit. We see thermalization occurring within microseconds. 



we will see the thermalization of both subspaces within 
microseconds. One can calculate the final occupation of 
the computational subspace from the initial density ma- 
trix (see description in Appendix A). For the time evo- 
lution of the excited state, we calculate Ooo « 0.82 and 
Poo « (0,0,0.44) (cf. Fig. 10(a)). A superposition of A 1 
and A'i evolves to « 0.70 and « (0, 0, 0.08) (cf. 

Fig. 10(b)). It is again seen explicitly that the long time 
evolution can be separated from the initial time evolution 
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(orange region). The initial time evolution on the inter- 
val [0, 10] ns can be described accurately by the second 
order Born approximation. 



VI. EFFECTIVE ERRORS 



We numerically simulate the time evolution of the 
qubit for different initial density matrices p (0) as in the 
previous section. We label the initial density matrix with 
the corresponding point P (0) = XV (crp (0)) on the Bloch 
sphere, with <r = (<7 x ,a y ,a z ). All analysis is done in the 
rotating frame with respect to % from Eq. (1). Again, 
we show the results from the numerical solution of the 
full master equation (cf. Eq. (56) with C rot — Cd from 
Eq. (16)) and compare them with the results of the cor- 
responding Nakajima-Zwanzig equation in second order 
Born approximation (Eq. (38) for the subspace qubit and 
Eq. (45) for the subsystem qubit). From these results we 
extract the error rates, as defined in Appendix B. We use 
the parameters of Sec. Ill A. The error rates extracted 
from the simulation of the full master equation and the 
second order Born approximation match to a high degree. 
We saw already in the analysis of the full time evolution 
in Sec. V, that on the interval [0, St] the time evolution 
of these equations match very closely. 

We will discuss the behavior separately for the initial 
points Pi = (0,0,1), P 2 = (1,0,0) and P 3 = (0,0,-1). 
We use in the following also the terminology "upper 
pole", "equator" and "lower pole" for these three points. 
We further need to extract the dephasing error on the 
equator. The trajectory starting at the upper and the 
lower pole will be restricted to the z-axis, which sets 

rf- P3 =0. 

To quantify the influence of noise on the qubit we need 
to compare the error rates 1^ 3 at Pj, as defined in Ap- 
pendix B, to the time of the experiment. The product 

yT^'StJ describes an error, which measures the leakage 

probability (t^ 3 Stj , the relaxation probability (v^ 3 8t\ 

and the dephasing probability (t^ 3 Stj . The entangle- 
ment fidelity F e from Eq. (B16) describes the effect of 
noise for all initial density matrices. We use the devi- 
ation of the entanglement fidelity from its ideal value, 
1 — F e , as a measure to quantify the overall error. 

We will find that all error rates can be described by 
only four toy models, introduced in Appendix C. For 
these models we can calculate the time evolution analyti- 
cally and describe with the extracted errors the expected 
error probabilities of the triple dot qubit. 




[A<eV] 



Figure 11. Energy diagram of the subspace qubit with major 
transition rates generated from local phase noise. The energy 
diagram shows only the relevant energy levels for this descrip- 
tion (cf. Fig. 2(a)). The transition rates can be grouped into 
three sets, which are described by the toy models analyzed in 
Appendix C. Model 1 describes pure relaxation of the qubit 
states (black arrows, see Appendix CI), model 2 pure de- 
phasing of the qubit states (gray arrows, see Appendix C2). 
Model 3 characterizes leakage of the qubit states to one state 
in the surroundings (orange arrows, see Appendix C3). 



A. Subspace Qubit 

1. Regime 1 

First, we analyze the subspace qubit in the (1,1,1)- 
regime (small e) with local phase noise from the inter- 
action with nuclear magnetic fields. As discussed in 
Sec. VA1, phase noise generates large error rates only 
for small energy differences (Regime 1 in Tab. III). Sim- 
ulations for local phase noise are shown in Fig. 12. We 
analyze phase noise on dot 1 and dot 2 separately and 
use Tf >2 = (20 ns) -1 (cf. Tab. II). 

The effective error rates can be understood when ana- 
lyzing the transition behavior in the DM. Quantum jump 
operations between the energy levels are only possible for 
the same s 2 -eigenstates. Since the interactions are local, 
the total spin quantum number S is not preserved. Ini- 
tially only the s z = \ doublet levels are occupied. Phase 
noise will mix within the s z = \ subspace (see sketch in 
Fig. 11). One should notice that all energy differences 
in the s z = \ subspace are small or comparable to the 
thermal energy. This makes transition rates for positive 
and negative energy differences similar. We can group 
all transition rates into three sets; each set corresponds 
to the error processes of a toy models from Appendix C. 
The transition rates of these toy models match to a high 
degree the results of the numerical solution of the full 
master equation, and of the calculation using the second 
order Born approximation. 

We want to discuss the results from Fig. 12 in detail 
and start with an analysis of local phase noise on the first 
quantum dot (cf. Fig. 12(a)). Model 1 of pure relaxation 
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Figure 12. Errors from local phase noise on the subspace qubit at E z = 100 /xeV. We take into account the influence of 
fluctuating hyperfine fields through the transition rates from Eq. (21). The effective error rates are extracted from the numerical 
simulation of the full master equation (blue lines) and of the second order Born approximation (green lines). Red lines describe 
effective errors from simple toy models, as described in the text. For each extracted error rate we add with a red number the 
specific model system under consideration from Appendix C. 



will mainly determine the Tj-behavior of the qubit (com- 
pare model 1 in Appendix CI). The interaction with 
hyperfine fields generates direct transitions between the 
qubit levels. These rates are large if the states are close 
in energy (compare h (er*,u;} in Eq. (21)). Hence the er- 
ror rates T^ 1 and T^ 3 vanish quickly when increasing the 



bias on the dots. Only a large difference of the transi- 
tion rate from the excited qubit level to the ground state 
compared to the reversed effect would cause large error 
rates T^ 2 and T^ 2 in model 1. Here both rates are either 
very similar (for small e) or both small (for finite e) . T^ 2 
and T^ 2 are therefore not described in the model of pure 
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relaxation. 

Transitions between the qubit subspace span { A 1 , A'i j 

and the quadruplet state Qi will determine leakage and 

also the error rate r p2 . We describe the effective error 
rates by the formulas from model 3 in Appendix C 3. We 
point out two important characteristics of this process. 
First of all, the transition rate from the A i -level to the 

2 

Qi -level is larger than the rate from the A 'i -level. The 

2 . 2 

smaller energy differences enhances the transition rates 
(compare h(al,uj) in Eq. (21)). Secondly, if 1 ' 2 and 
2 is larger for positive bias than for negative one. This 



results because the transition amplitude j^j A: 

is larger at positive bias than at negative bias. The tran- 
sition amplitude can be read off directly when comparing 
the eigenstates for e = ±oo with the eigenstate at e = 
in Fig. 7. 

The error rate r p2 has additionally a very interesting 
behavior. It can be described mainly by pure dephas- 
ing (see model 2 in Appendix C2). We especially find a 
point where the dephasing rate has a minimum. When 
the energy level fluctuations at both levels are equal, r p2 
vanishes (cf. Eq. (Cll)). We can determine this point 

analytically ^Ai \a\\ Aij = ^A'i \a\\ A'^ and find 

J12 = 2J 2 3. For the chosen parameters in our calcula- 
tion, we can approximate e ~ — e ^-. 

For phase noise on dot 1 the behavior of r p3 is not cap- 
tured by the error rates from all three error models. This 
is caused by the large dependency of the transition rates 
on the energy differences (compare h {a\, to) in Eq. (21)). 

Instead of a direct transition ^A'i — > Q^j, we observe 

rather a two step process ^A'i — > Ai — > Qi . 

We can summarize the contribution of all error rates 
to the deviation of the entanglement fidelity from its 
ideal value (1 — F e ). The overall behavior is determined 
mainly by the dominant error rate, which is in this case 
dephasing at the equator. We will however never find a 
value e for which the errors due to phase noise on dot 1 
have a vanishing effect. 

Local phase noise on dot 2 does not allow any transi- 
tions from the A i -eigenstate if there is no external bias. 
This can be understood when analyzing the eigenstates 
of Hamiltonian (1) at e = in Fig. 7. The Ai-state in- 
volves a singlet state on the outer dots, while the remain- 
ing eigenstates contain only triplet states. Since interac- 
tions on the middle dot leave the states on the two outer 
dots untouched, the A i -state is protected from any local 
noise on the middle dot. This directly forbids leakage 
and relaxations from the upper pole at e = 0. Further, 

3 (e = 0) = 0, since the lower pole never relaxes to the 
upper one. 

The remaining features of the transition rates can be 
understood from their strong energy dependence (com- 
pare h{a 2 ,u) in Eq. (21)). Transitions from the qubit 
states to the Qi -quadruplet states describe leakage er- 



rors r^' P2 ' P3 through model 3 (cf. Eq. (C13))-(C15)). 
Leakage from the A 'i -state (r Pa ^ decreases strongly at 
finite bias, since the energy difference to the Qi -state in- 
creases. The leakage error from Ai (r Pl ^) has the oppo- 
site characteristic. T p2 and T p2 are mainly determined 
by the average value of T Pl and T p3 . Since T Pl » T p3 , 

we get approximately r o 2 = T-^ 2 « — ^— . 

Model 1 describes relaxation of the states at the upper 

and the lower pole ^r Pl,P3 ^. The energy difference of 

the qubit states Ai and A'i is minimal at e = (cf. 

2 2 

Fig. 11). The relaxation rate h (cr^w) is large for small 
e. Only the special symmetry at e — causes the large 
decrease of r Pl,r>3 directly around e = 0, where the tran- 
sition amplitude vanishes. 

Dephasing errors ^r p2 <5i^ are describing by the pure de- 
phasing mechanism of model 2 (see Appendix C2). We 
can calculate the difference in the energy level fluctuation 
of the two states: 



Ai 



A, 



A\ \a 2 \ A\ 



3 \J Jf 2 ~ Ji2-hz 



T 2 



(57) 



Eq. (57) has the limit | for Ji2< J23 and | for J 12 = J 2 3- 
These limits determine for the most part the overall error 
of the qubit (1-F e ). 

We find that local phase noise induces large errors to 
the time evolution of the subspace qubit. Especially pure 
dephasing, as described by toy model 2 (see Appendix 
C2), limits the performance of the triple quantum dot. 
Large errors are generated via phase noise on dot 1 for 
strong external bias, while phase noise on dot 2 is critical 
for the unbiased dot. This effect can be understood when 
considering the high symmetry regimes of the Hamilto- 
nian in Fig. 7. Phase noise is always most critical when 
it acts on the eigenstates of a single quantum dot. 



2. Regime 2 

Next we analyze errors of the subspace qubit though 
local spin relaxation in the low bias regime (e small) . We 
consider only electric bias e in the range [e_,e+]. Lo- 
cal spin flip operators can generate transitions, changing 
the angular momentum quantum number (As z = ±1). 
The total spin quantum number S is not necessarily pre- 
served. Transition rates through hyperfine interactions 
are highly suppressed due to the large energy difference 
of states with different s^-quantum number (cf. Fig. 2(a) 
and the transition rates in Eq. (21); the main effects are 
captured in Regime 2 of Tab. Ill) . For phonon- mediated 
transitions the level splitting must be quite large to see 
transitions in the nanosecond regime (cf. Eq. (27)). At 
E z = 100 /j,eV we only see the effect of two transitions 
rates. All other transition rates are greatly suppressed 
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Figure 13. Description of local spin relaxation for the sub- 
space qubit at large external magnetic fields (E z = 100 /xeV). 
Only local spin relaxations from the interaction with phonons 
are significant, (a) Sketch of energy diagram of the subspace 
qubit to describe the time evolution under local spin relax- 
ation. Ti and To are the dominant transition rates from the 
qubit levels to the s z = | quadruplet state, (b) Transition 
amplitudes for the qubit state i through the noise operator 
<t + acting on dot j: (crj_) . = j 3 \<j j + \ Wij\ (Wi — Ai, 
Wo = A'x). 



(cf. Fig. 13(a)). In Fig. 14 we show the error probabili- 
ties for this parameter regime. 

To explain the error probabilities in detail, we only 
need to discuss two properties of the relevant transition 
rates. First of all we need to analyze the transition am- 
plitude from the qubit levels |Ai, A'i j to the s z = |- 

quadruplet state. These amplitudes are drawn as a func- 
tion of the external bias e in Fig. 13(b). For relaxations 

3 I oil A] 



on dot 1 the transition amplitude 

creases steadily from negative to positive bias. For the 

relaxations from the A \ -state this effect is reversed. The 

2 

transition amplitudes are equal at J12 = 2J23. One can 



proof this immediately, when looking at the eigenstates 

Ai and A'i of the triple dot Hamiltonian H in Eq. (7)- 

2 2 

(8). For spin relaxation on dot 2 the transition ampli- 
tude from A'i is always greater than the transition am- 

2 

plitude from A 1 . Additionally, we see a maximum of 



Qa I oil 



Ai 



vanishes 



at zero bias. 

at this point. Second, the transition rates depend on 
the energy difference between the quadruplet state and 
the doublet states. We won't see any effect from transi- 
tions involving the A 1 -state, since its eigenenergy is only 

weakly dependent on e. For the A'i -state the eigenenergy 

2 

is influenced strongly by e. 

Having these two discussions in mind we can under- 
stand the effective error rates of Fig. 14. We only need 
to compare the results with the analysis from model 3 in 
Appendix C 3. Here the special case of transitions from 
the qubit to the surroundings, without its opposite effect, 
applies (cf. Eq. (C20)-(C26)). 

Leakage at the upper pole (r^ 1 ^ is only dependent 
on the transition rate from Ai to Qa (cf. Eq. (C20)). 
For relaxation at dot 1 the leakage rate therefore steadily 
increases, while relaxation on dot 2 has a local minimum 
at zero bias. For the lower pole P3 the effect is reversed. 
Leakage on dot 1 steadily decreases with e, while it has 
a maximum at e = for noise acting on dot 2. We 
also see the dependence on the transition energy at large 
bias. Here the error rate T^ 3 decreases further. The 

leakage rate at the equator (r^ 2 ^ is represented in lead- 
ing order by the average leakage rates at Pi and P3 (cf. 
Eq. (C21)). 

Relaxation at Pi and P 3 vanishes in model 3 in lead- 
ing order (cf. Eq. (C23) and (C25)). The results from 
the numerical analysis are not obtained in model 3. But 
note that their magnitude is very small. Relaxation at 

the equator \ Ff is determined in leading order by the 

difference in transition rates from A 1 to Q 3 , compared 

to the rate from A\ (cf. Eq. (C24)). We see vanish- 

2 

ing relaxations at Jyi — 2J23 for noise acting on dot 1. 
Dephasing at the equator T^" 2 shows a very similar char- 
acteristic. It is also dependent on the difference in the 
transition rates (cf. Eq. (C26)). 

The detected entanglement fidelity for local spin relax- 
ation is close to 1. For local noise on dot 1 or 2 it does not 
show a characteristic dependence on the bias parameter. 
In total all resulting errors are much smaller compared to 
the influence of phase noise discussed in the last section. 



B. Subsystem Qubit 

The subsystem qubit is, for large external magnetic 
fields, equivalent to the subspace qubit. We won't add 
any results for the subsystem qubit in this case. As de- 
scribed in Sec. IV B 3, phase noise also acts on the sub- 
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Figure 14. Errors from local spin relaxation generated from interactions with phonons for the subspace qubit in the (1, 1, 1)- 
regime at E z — 100 /xeV. Blue lines are calculated from the numerical simulation of the full master equation. Green lines are 
obtained from the second order Born approximation. The red lines represent the results from the analysis of model 3, which 
involve only two transition rates (compare description in the main text). We can see that all descriptions are matching well. 



system qubit identically at small external magnetic fields. 
The subsystem qubit differs from the subspace qubit only 
at small external fields and for local spin relaxations. 
We analyze the subsystem qubit only for magnetic field 
strengths comparable to the thermal energy. Here both 
the s z = \ and s z = —^-subspace is initially occupied 
(see distribution function in Eq. (12)). The qubit levels 
are characterized by the formally introduced ^-quantum 
number, as described in Sec. II B. 



1. Regime 3 

Already in the small detuning regime, levels of dif- 
ferent s z -quantum number cross. Hyperfine interactions 
can generate transitions between these levels through lo- 
cal spin flips. The total spin quantum number S is not 
preserved by local interactions. In Fig. 15 we extract the 
error rates generated by local spin relaxation on dot 1 
and 2. Relaxation due to phonons won't be detectable 
for the subsystem qubit in the regime of nanoseconds. 
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Larger energy differences are needed to see strong effects 

(see transition rates in Eq. (27)). 

Phase noise from hyperfine interaction generates large 

transition rates only if energy levels of different s z - 

quantum number are close in energy (compare h (<r* , oS) 

in Eq. (21); Regime 3 in Tab. III). This is the case at 

the level crossings of the energy diagrams (cf. Fig. 16). 

Essentially two sets of transitions are important. First 

of all, leakage of the doublet levels to quadruplet states 

plays an important role (model 3 from Appendix C3). 

Since local spin flip operators always change the s z - 

quantum number, only transitions from the I = 1 states 

to the quadruplet levels are possible. Leakage from the 

I = states is highly suppressed. Their energy differs 

significantly from the quadruplet levels of different s z - 

quantum number. Second, there are internal transitions 

between the subsystem states Ai and A'_ i , as described 

2 2 
with model 4 from Appendix C 3. 

All effects can be summarized easily. Model 4 deter- 
mines the relaxation errors at the upper and lower pole 
as well as the dephasing properties. Model 3 describes 
the leakage behavior. Since the internal transition rates 
are very similar at the equator, effectively no relaxation 
is generated internally by transitions between the Ai 
and A'_i states. Relaxation is rather determined by the 

2 

indirect process of leakage to the quadruplet levels. 

For relaxations acting on dot 2 we see the symmetric 
error rates at positive and negative bias. For relaxations 
on dot 1 the error rates for positive bias is greatly sup- 
pressed. A'i approaches for positive bias the eigenstate 

2 

from e — > oo in Fig. 7. It involves a singlet state on the 
strongly coupled dots 2 and 3. Noise on dot 1 leaves 
this singlet untouched. This state does not couple to any 
quadruplet state or I = state. 

Finally we want to explain why Fj^ 1 is not correctly 
described by our simple model analysis. Leakage at the 
upper pole is determined by a second order process of 
internal transitions followed by leakage (Ai ^ A'_i — > 

Q\ ) . We describe this process neither by model 3 nor by 

2 

model 4. 

Overall we find that at the point of level crossings fluc- 
tuating hyperfine fields can generate major errors of the 
subsystem qubit. The asymmetry of the error rates for 
spin relaxation on dot 1 between positive and negative 
bias is a very interesting result. For positive bias fluctu- 
ating hyperfine fields generate nearly no errors, since the 
Hamiltonian (1) is approaching a high symmetry regime 
as discussed in Sec. IV B 4. 



VII. SUMMARY AND OUTLOOK 

The exchange-only qubit has been analyzed with two 
different schemes for implementing a coded qubit in the 
Hilbert space of three single occupied quantum dots, re- 
sulting in the subspace and the subsystem qubit. The re- 
laxation and decoherence dynamics of both these qubits 



has been calculated, with both nuclear spins and phonons 
taken into account. These interactions are described in 
the Davies model (DM), a particular Markov approxi- 
mation with a transparent quantum-jump interpretation 
and consistent long-time behavior. The systematics of 
the early-time dynamics, which is of most interest for 
qubit experiments, is quite distinct from the systemat- 
ics of the long time evolution. We have focused on the 
initial time evolution and have extracted errors for the 
subspace and the subsystem qubit. We can describe all 
results by relating them to just four toy models, whose 
time evolution can be calculated analytically. 

For local phase noise, arising from the interaction with 
fluctuating hyperfine fields, the influence on the the sub- 
system and the subspace qubit are identical. Local phase 
noise is critical for GaAs systems; it is the strongest 
mechanism for the loss of phase coherence. Sizable phase 
errors for both the subspace and subsystem arise after 
just 10 nanoseconds of evolution. 

The influence of local spin relaxation is different for the 
subspace and the subsystem qubit. Since the subspace 
qubit is always operated at large external magnetic fields, 
spin relaxations from the interactions with phonons need 
to be considered. This effect generates large transition 
rates only between energy levels with large energy differ- 
ence. Our analysis shows that in GaAs systems, oper- 
ated at large external magnetic fields, only small errors 
are generated. For the subsystem qubit an operation at 
small magnetic fields is also possible. Here, only local 
spin relaxations from the interactions with nuclear mag- 
netic fields are important. These interactions generate 
large errors at the crossing of energy levels of different 
s z -quantum number. This process has a very interest- 
ing property for phase noise acting locally on one of the 
outer quantum dots. Errors can be highly suppressed 
depending on the sign of the bias parameter e. 

To state our results briefly, our analysis shows that in 
GaAs samples (large nuclear bath) the subsystem and the 
subspace qubit have about equally good coherence prop- 
erties. Both qubit implementation schemes suffer from 
local phase noise, generated from fluctuating hyperfine 
fields. Spin relaxation from hyperfine fields will be im- 
portant only at the point of level crossings. If these points 
can be avoided when manipulating the qubit, spin relax- 
ations induced by fluctuating hyperfine fields are negli- 
gible. If one attempts to use the crossing points in the 
energy level diagram for qubit manipulations (cf. the 
attempt to manipulate a singlet-triplet qubit at cross- 
ing points in the energy diagram 37 ), one has to pay at- 
tention to fluctuating hyperfine fields. Interactions with 
phonons will usually be less important. This mechanism 
will only be significant if there are strong external mag- 
netic fields and large energy differences. The interaction 
with phonons can completely depopulate the qubit, but 
in GaAs systems this evolution only occurs on the mi- 
crosecond timescale. 

Overall it is a very interesting task to test the local 
nature of the error models. Especially for the influence 
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Figure 15. Errors of the subsystem qubit generated by local spin relaxations for small external magnetic fields (E z = 2.5 /ieV). 
Large error rates are only seen at the point of level crossings. For noise acting on dot 1 the errors at positive detuning are 
highly suppressed. Blue lines represent the calculation via the simulation of the full master equation. Green lines represent the 
second order Born approximation. Red lines show results from the analysis of model system 3 (see Appendix C 3) and model 
system 4 (see Appendix C4). 



of nuclear spins, which behave on short time scales like 
classical fluctuating magnetic fields, the local influence 
of the qubit dynamics is worth testing. Such an experi- 
ment would require the control of the randomness of the 
hyperfine fields at the positions of the different quantum 
dots. If it is possible to reduce the randomness at two 
of the three quantum dots, so that the hyperfine interac- 
tion noise acts dominantly on one of the three quantum 
dots, one can try to test the different scaling behavior of 



the error rates with the bias parameter e. Furthermore 
our analysis method in the DM should be helpful for the 
description of other coded qubits implemented in more 
complex Hilbert spaces. We show in detail here, for the 
triple-dot qubit, how the interaction with complicated 
baths can be reduced to just an effective evolution on 
the coded qubit itself. Such analysis could be extended 
to other coding strategies, when the need arises. 
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Appendix A: Long Time Limit of Time Evolution 
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Figure 16. Transition diagram for subsystem qubit when 
hyperfine interactions generate local spin flips. Large error 
rates are observed only at the region of level crossings. They 
can be described by leakage transitions to a quadruplet state 
(model 3 from Appendix C 3) or internal transitions between 
two states of the subsystem qubit (model 4 from Appendix 
C4). 



Since in the DM the system equilibrates to thermal 
equilibrium, we can calculate the long time behavior of 
the models of Sec. V analytically. One needs to pay atten- 
tion that only subspaces that are connected by internal 
transitions equilibrate. 



For the subspace qubit under the influence of phase 
noise (see analysis in Sec. V A 1) , we can restrict ourselves 

to the subspace |q i , A i , A'i j. In the long time limit, 

the density matrix will show partial equilibration: 



, I 2 2 2 J 



1/2 



/ E 

ie { Q §' A i' A *} 



(Al) 



The long time limit for the population of the subspace 
qubit can be obtained from the total leakage to the 
quadruplet state: 



B 



^1/2 



(A2) 



1/2 



Since the off-diagonal elements of the density matrix van- 
ish, it is clear that = 0. The long time limit of the 
qubit 's polarization can be calculated from the difference 
in the population of the Ai and A 'i -state: 



'-Q1/2 



(A3) 



n 2 



+ e 



Eq. (A2) and (A3) are used to calculate the long time 
limit for the qubit evolution in Sec. V A 1 (see especially 
insets of Fig. 8). 

The subsystem qubit with local spin relaxations has 
a description which is slightly more complicated. The 
simulation of Sec. V B 1 analyzes the specific situation of 
phase noise near the crossing points of energy levels (see 



orange line in energy diagram Fig. 2(b)). We take into 
account only transitions that occur on the time scale of 
microseconds. This limits us to transitions in two sub- 
spaces sspl and ssp2: 



.-pl = {0j,Aj,A^} 
ssp2 = {q§,A'i} 



(A4) 
(A5) 



In both subspaces thermal equilibrium is reached. Tran- 
sition rates between these two subspaces and to the re- 
maining states are very small. Only the I = and I = 1 
states are occupied at t — for a subsystem qubit (see 
qubit definition in Sec. II B). The final values for the 
qubit evolution are dependent on the initial density ma- 
trix (cf. Eq. (12)): 



P(°)= p p )®P° 



(A6) 



p (0) determines the part of the density matrix, which is 
initially part of sspl (O sspl ), of ssp2 (O ssp2 ) or remains 
unchanged (O u = 1 - O sspl - O ssp2 ) . The initial pop- 
ulation of subspace sspl depends on the occupation of 
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the states Ai and A' t . It can be described by the en- 

2 ~2 

tries Pn and P 00 of p(0) from Eq. (A6), which are itself 
related to the initial polarization P z (0): 



QSSpl 



1 + P z (0) 1 
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z 1 + e T * 
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z 1 + e t k 
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For subspace ssp2 only the initial occupation of the state 
A'i plays a role, which leads to 



0ssp2 



1 - P z (0) 



1 



1 + e T *r 



(A8) 



The final population of the subsystem qubit is deter- 
mined by all transition rates to the quadruplet states: 
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All superpositions vanish in the long time limit (X^ = 
0) and the final polarization can be calculated from the 
difference in population of the I — states and the I = 1 
states: 



Zqo — O; = o — Oi—i 
=O u + O sspl - 
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Eq. (A9) and Eq. (A10), together with Eq. (A7)-(A8), 
can be used to describe the long time limit for the subsys- 
tem qubit in Sec. VB 1 (see especially insets of Fig. 10). 



Appendix B: Analysis Model for Single Qubit 
Description 

In this section we connect the description for the sin- 
gle qubit time evolution from the spin-based quantum 
computation community to the common one in quantum 
information theory. In the first one, the qubit evolution is 
described by the evolution on the Bloch sphere (compare 
e.g. the recent review of Kloeffel and Loss 38 ). One com- 
monly uses maps on density matrices in an information 
theoretical approach 39 . 



1. Solid State Approach 

In a solid state approach, one commonly uses two spe- 
cific time scales to describe the evolution on the Bloch 
sphere, which originally came up in the literature of 
NMR 40 . First the longitudinal relaxation time T\ de- 
scribes the evolution from the excited qubit state 1 1) to 
the ground state |0). We will call this time scale "re- 
laxation time" in the following. Second the transverse 
relaxation time T 2 (which we will call "dephasing time") 
describes the relaxation of a quantum mechanical super- 
position (| 1) + 1 0) ) /y/2 to a mixed state. 

We describe a complex time evolution in our analy- 
sis, including leakage from the computational subspace 
to the embedding Hilbert space. We characterize this 
evolution by the introduction of a third time scale, which 
we call "leakage time" T . Even though all parameters 
are originally meant to describe the inverse rates of expo- 
nential time evolutions, we are fitting our results of more 
complex dynamics to these parameters. We analyze the 
initial time evolution from points P (0) = Tr (erp (0)) 
on the Bloch sphere, with er = (a X7 a y7 a z ), and extract 

the leakage rate r^ >( ' \ the relaxation rate T^ -* and the 

dephasing rate T^ ^ of the initial time evolution. We 
correct all rates by a factor linear in the time argument 
to help account for the non-exponential behavior: 



^P(O) P(0) . P(0)jf. 
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The leakage time T 



P(0) 



is described by the corresponding leakage rate r o 



P(0) 



(A10) (Tq( ^ of the trace evolution of the relevant part of 



the density matrix (see description in Sec. IV A): 



Pp(o)(«)«(i-if (0) «*+(r? (0) ) 



2^ 
2 



(B2) 



Since leakage leads to a depopulation of the qubit, we 
renormalize all Bloch sphere parameters by the trace 
of the relevant part of the density matrix (Pi (t) = 
Pi(t) /Tr (p rel (t))). We assign the relaxation time 
j,P(0) _ ^pP(o)^ ^ ^ e z-evolution of the qubit from 
the initial polarization P z (0) to the final polarization 



p z (5t) «p z (o)(i-rf (0) 5t+(r[ 



-.P(O) 



2 St 2 



(B3) 
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Dephasing describes the loss of phase coherence of a 
qubit. We especially refer to the relaxation of quan- 
tum mechanical superpositions to a mixed state. On the 



25 



Pi = (0,0 



Pi 
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Figure 17. Sketch of transition rates, extracted from the 
initial time evolution of the qubit on the Bloch sphere at three 
special points. The Bloch sphere parameters are renormalized 
by the trace evolution of relevant part of the density matrix, 
which describes the qubit: P XiZ (t) = P XlZ (t) /Tr (p rel ). 



Bloch sphere it is connected to the rate at which a point 
on the surface of the Bloch sphere relaxes to the z-axis. 



We extract the dephasing time T 2 
from the initial time evolution: 



P(0) 



2T 



P(0) 



41 



actly follows the z-axis (cf. Sec. IV B 2). From one point 
of the equator, e.g. P (0) = P2 = (1, 0, 0), we extract the 
leakage, relaxation and dephasing rates. The relaxation 
rate is extracted from the time evolution to the north or 
the south pole. Initially just one of the rates, defined in 
Eq. (B4), is positive. This positive number defines the 



relaxation rate r. 



A sketch of all transition rates on 



the Bloch sphere is shown in Fig. 17. 



2. Information Theoretical Approach 

Commonly, one describes the time evolution through 
the completely positive linear map est in an information 
theoretical approach, est constructs from the initial den- 
sity matrix p(0), the density matrix at some later time 
p{5t): 



p(St)=est (p(0)) 



(B5) 



In our analysis est is trace decreasing, since we also 
take into account leakage to the surroundings. Due to the 
special trajectory generated in our model (cf. Sec. IV B), 
we only need 7 free parameters to completely describe 
the initial time evolution of our system. Since the map 
est is linear, it is fixed completely by its action on four 



pure states 1 0) , 1 1) and 



1) 



1 

(i) 



0) /V2: 



P x (St) 



i-r 



P(0) 



5t 



,P(0) 



2^2 
~2~ 



P x (0) (B4) 



Following the discussion of Sec. IV B 2, we can restrict 
the analysis to just one plane (e.g. x-z-plane). It is suf- 
ficient to extract the parameters r^' ', T^ ^ and T^ ^ 
only for three values of P (0), i.e. only at three points on 
the surface of the Bloch sphere, to describe the full time 
evolution of the qubit (cf. Sec. B 2). From the upper and 
the lower pole P (0) = Pi(3) = (0,0, (-)l), we extract 
the leakage and the relaxation rate to the opposite pole. 
Because of the properties of the DM, the trajectory ex- 



e (|l)(l|)=a 1 |l)(l|+a 2 |0) (0| 
e(|0) (0|)=O3 |1)(1| +a 4 |0) (0| 

=a 5 |1)(1| +a 6 |0) (0| 



a 7 



(B6) 
(B7) 

(B8) 



It is straightforward to relate the parameters ai — 07 
to the evolution rates T t 3 , that were defined earlier (cf. 
Eq. (B2)-(B4)): 
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Various properties of interest can be calculated from system R: p R ® = |m) /2. It describes how well 

the map est- One example is the entanglement fidelity: the entanglement between two systems is preserved under 

F e = Tr (p RG (1 ® e$ t ) (p^)). is the maximally the action of the noisy quantum channel est 42 - 

entangled state of the noisy system Q and the reference 
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Appendix C: Model Systems 



r 



toy model systems individually. 



We present in Fig. 18 four model systems to describe 
the effective error rates of the qubit defined in triple 
quantum dot systems. In the DM energy eigenstates can 
couple through quantum jumps. Transitions are possi- 
ble between the qubit levels 1 1) and 1 0) , but also to 
other states of the embedding Hilbert space | Out) . The 
strength of the quantum jumps is specified by the coef- 
ficients Tj e M, that are determined by two constants. 
First of all, there are the transition rates h (A, uj), which 
are extracted from experiments (for the interaction with 
hyperfine fields, cf. h (a l z /a x , uS) in Eq. (21); for the in- 
teraction with phonons, cf. h{a l x ,uj) in Eq. (27)). Sec- 
ondly, we need also the matrix elements of the transition 
operator between the energy eigenstates (cf. Eq. (18)). 
While Tj can be negative, only the positive number T 2 
describes a rate. In the following, we describe all four 



1. Model 1: Pure Relaxation 

The model of pure relaxation describes the transitions 
in a two level system through raising and lowering oper- 
ators: T ± er ± (cf. Fig. 18(a)). The full master equation 
is: 

P (t)=V [T + a + ] (p (i)) + V [T_a_] (p (i)) (CI) 
= T 2 + V[a + }(p(t)) + T 2 _V[a.](p(t)) (C2) 

Eq. (C2) can easily be solved. The effective error rates, 
defined in Sec. B, can be extracted. Assuming T 2 _ > T 2 + , 
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(a) Model 1: Pure Relaxation (b) Model 2: Pure Dephasing 




It is clear that no leakage arises in this model. Relax- 
ation from the Pi and P3 is determined by the direct 
transitions to the opposite pole. The combination of the 
two transition rates gives a reduction of the overall error. 
At P2, we initially see relaxation to the lower pole with 
a rate determined by the difference of the two transition 
rates. 



2. Model 2: Pure Dephasing 
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(c) Model 3: Two State Leakage (d) Model 4: Internal 

Transitions of the Subsystem 
Qubit 

Figure 18. The time evolution of the qubit in triple quantum 
dots can be described by four toy models (a)-(d). 



Pure dephasing (cf. Fig. 18(b)) is described in the DM 
by the coupling operator: 



A = 




(C8) 



When solving the master equation p(t) = T> [A] (p(t)), 
we extract the following transition rates: 
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The coupling operator (C8) generates neither relaxation 
nor leakage. A describes fluctuating energy levels, which 
leads to pure phase noise. 



we g e |.. 3. Model 3: Two State Leakage 

In our calculation, we need to describe leakage of the 
u v^°) qubit states to exactly one state of the embedding Hilbert 

space (cf. Fig. 18(c)). When solving the master equation 

p(t) =V [T 1+f j^° ut ] (p(i)) +2? [T^a^*] (p(t)) 
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we can extract the effective error rates (assuming T 2 , > 
T u+ ): 
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We need in particular two special cases of this model. 
First of all, for the subspace qubit we analyze the case 
where only the transition rates from the qubit levels to 
the surroundings are significant. We can set T 1+ = T 1; 
T + = T 2 and Ti_ = Y - = and obtain the error 
rates: 
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On the north pole Pi of the Bloch sphere neither leakage 
nor relaxation is seen. On the south pole P3, we observe 
pure leakage. 



Model 4: Internal Transitions of the Subsystem 
Qubit 



Here, for the north and south pole (Pi and P 3 ) no re- 
laxation is generated and only leakage occurs. 

Second, for the analysis of the subsystem qubit we 
need to describe the transition of only one qubit level 
to the surroundings. For T 0+ = T + and T _ = T_, 



For the analysis of the subsystem qubit, we will need 
the extract error rates at the crossing of two levels defin- 
ing the qubit. We can solve the Davies master equa- 
tion p(t) = V[A](p(t)), describing the toy model of 

Fig. 18(d). For the case T 2 > e"^T^, we get the 
error rates: 



29 



r Pl,P 2 ,P3 
1 n 



t2 
1 10 



t2 
1 10 



Ez 



-T?„ 



~v~2 
1 10 



r ^ 1 



e 1 + e 

2 



1 



' -'Til! + / Tf -e-#Y^ ' 



1 + e t k 



1 



e 

l + e~^ 1 

"2 i — t^t2 



*Y^2 

5^ 1 01 



T§1 



ft 



p 2 _Tf + e Tqj 



it 2 = 



2 1 



(Tf - Tgi) 



, v^2 

3^ 1 01 



2 ft 



ft 



f 01 



(C34) 
(C35) 

(C36) 

(C37) 

(C38) 



s.mehl@fz-juelich.de 

d .divincenzo@fz-juelich . de 

R. Hanson, L. P. Kouwenhoven, J. R. Petta, S. Tarucha, 
and L. M. K. Vandersypen, Reviews of Modern Physics 
79, 1217 (2007). 

D. Loss and D. P. DiVincenzo, Physical Review A 57, 120 
(1998). 

J. Levy, Physical Review Letters 89, 147902 (2002). 
J. M. Taylor, H.-A. Engel, W. Diir, A. Yacoby, C. M. Mar- 
cus, P. Zoller, and M. D. Lukin, Nature Physics 1, 177 
(2005). 

J. R. Petta, A. C. Johnson, J. M. Taylor, E. A. Laird, 
A. Yacoby, M. D. Lukin, C. M. Marcus, M. P. Hanson, 
and A. C. Gossard, Science 309, 2180 (2005). 

A. C. Johnson, J. R. Petta, J. M. Taylor, A. Yacoby, M. D. 
Lukin, C. M. Marcus, M. P. Hanson, and A. C. Gossard, 
Nature 435, 925 (2005). 

D. P. DiVincenzo, D. Bacon, J. Kempe, G. Burkard, and 
K. B. Whaley, Nature 408, 339 (2000). 

E. A. Laird, J. M. Taylor, D. P. DiVincenzo, C. M. Marcus, 
M. P. Hanson, and A. C. Gossard, Physical Review B 82, 
075403 (2010). 

L. Gaudreau, G. Granger, A. Kam, G. C. Aers, S. A. Stu- 
denikin, P. Zawadzki, M. Pioro-Ladriere, Z. R. Wasilewski, 
and A. S. Sachrajda, Nature Physics 8, 54 (2012). 
L. Viola, E. Knill, and R. Laflamme, Journal of Physics 
A 34, 7067 (2001). 

B. H. Fong and S. M. Wandzura, Quantum Info. Comput. 
11, 1003 (2011). 

J. R. West and B. H. Fong, arXiv: 1203.4296 [quant-ph] 
(2012). 

J. Kempe, D. Bacon, D. A. Lidar, and K. B. Whaley, 
Physical Review A 63, 042307 (2001). 

D. A. Lidar and K. B. Whaley, Lecture Notes in Physics 
622, 83 (2003). 

E. B. Davies, Communications in Mathematical Physics 
39 , 91 (1974). 

E. Fick and G. Sauermann, Quantenstatistik dynamischer 
Prozesse (Harri Deutsch, 1983). 

E. Fick and G. Sauermann, The Quantum Statistics of Dy- 
namic Processes (Springer, 1990). 



18 W. A. Coish and D. Loss, Physical Review B 72, 125337 
(2005). 

19 J. M. Taylor, J. R. Petta, A. C. Johnson, A. Yacoby, C. M. 
Marcus, and M. D. Lukin, Physical Review B 76, 035315 
(2007). 

20 These parameters are related to the ones from the Hubbard 
Hamiltonian in Eq. (2) by e_ = — (ei — 62 + Ui) and e+ = 
£3 — £2 + Us. 

21 M. Celio and D. Loss, Physica A 158, 769 (1989). 

22 S. Bravyi and J. Haah, arXiv:1112. 3252 [quant-ph] (2011). 

23 H. Spohn, Letters in Mathematical Physics 2, 33 (1977). 

24 We refer especially to typical qubit manipulation times. In 
the mentioned publication by Gaudreau et al. pulse times 
well below 10 ns are used 9 . 

25 In the paper by Gaudreau et al. the tunnel couplings T are 
defined differently than the parameters t in Eq. (3) and (4) . 
The two constants are however connected by T — t/y/2. 

26 W. A. Coish and D. Loss, Physical Review B 70, 195340 
(2004). 

27 W. A. Coish, J. Fisher, and D. Loss, Physical Review B 
81, 165315 (2010). 

28 L. Cywinski, W. M. Witzel, and S. D. Sarma, Physical 
Review Letters 102, 057601 (2009). 

29 L. Cywinski, W. M. Witzel, and S. D. Sarma, Physical 
Review B 79, 245314 (2009). 

30 T. D. Ladd, Physical Review B 86, 125408 (2012). 

31 I. A. Merkulov, A. L. Efros, and M. Rosen, Physical Re- 
view B 65, 205309 (2002). 

32 In single spin experiments fluctuating in-plane magnetic 
fields, both in x- and y-direction, cause spin flips. The gen- 
erated dynamics is very similar. We will just consider the 
coupling operators a l x , since they directly relate to single 
spin flips a±. 

33 S. Amasha, K. MacLean, I. P. Radu, D. M. Zumbiihl, M. A. 
Kastner, M. P. Hanson, and A. C. Gossard, Physical Re- 
view Letters 100, 046803 (2008). 

34 K. C. Nowack, M. Shafiei, M. Laforest, G. E. D. K. 
Prawiroatmodjo, L. R. Schreiber, C. Reichl, W. Wegschei- 
der, and L. M. K. Vandersypen, Science 333, 1269 (2011). 

35 X. Hu and S. D. Sarma, Physical Review Letters 96, 
100501 (2006). 



30 



K. Blum, Density matrix theory and applications (Plenum, 
New York, 1996). 

H. Ribeiro, J. R. Petta, and G. Burkard, Physical Review 
B 82, 115445 (2010). 

C. Kloeffel and D. Loss, arXiv:1204.5917 [cond-mat.mes- 
hall] (2012). 

M. A. Nielsen and I. L. Chuang, Quantum computation and 
quantum information (Cambridge University Press, Cam- 



bridge, 2000). 

C. P. Slichter, Principles of magnetic resonance (Springer, 
Berlin, 1990). 

The factor 2 is used to agree with the definitions of Nielsen 
and Chuang 39 . 

D. C. Marinescu and G. M. Marinescu, Classical and quan- 
tum information (Elsevier Academic Press, Amsterdam, 
2012). 



